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Topological behavior can be masked when disorder is present. A topological insulator, either 
intrinsic or interaction induced, may turn gapless when sufficiently disordered. Nevertheless, the 
metallic phase that emerges once a topological gap closes retains several topological characteristics. 
By considering the self-consistent disorder-averaged Green function of a topological insulator, we 
derive the condition for gaplessness. We show that the edge states survive in the gapless phase 
as edge resonances and that, similar to a doped topological insulator, the disordered topological 
metal also has a finite, but non-quantized topological index. We then consider the disordered 
Mott topological insulator. We show that within mean-field theory, the disordered Mott topological 
insulator admits a phase where the symmetry-breaking order parameter remains non-zero but the 
gap is closed, in complete analogy to 'gapless superconductivity' due to magnetic disorder. 



I. INTRODUCTION 

The discovery of topological insulators (TI) has 
brought us to nothing less than a revision of electronic 
band structure theory. Topological behavior was pre- 
dicted for two-dimensional (2d) systems such as graphene 
with spin-orbit coupling 1 ' 2 and HgTe quantum wells 3 . 
While the spin-orbit coupling in graphene turns out to 
be too weak, topological behavior has indeed been ob- 
served in HgTe/CdTe heterostructures 4 . Shortly there- 
after, three-dimensional (3d) topological insulators were 
predicted to exist in bismuth alloys 5-7 , and by now have 
been discovered in a variety of materials (see, e.g., Rcfs. 
8 and 9). From a theoretical perspective, these mate- 
rials are described by a non-zero quantized topological 
index, which is essentially a generalization (albeit a non- 
trivial one) of the Chern number 5-7 ' 10-12 . As in Hall 
insulators, the topological index implies the existence of 
protected edge/surface states at the boundaries of the 
system 10 ' 13 , which are responsible for a variety of excit- 
ing aspects of topological insulators. In two dimensions, 
the Id edge states yield a quantized two-terminal conduc- 
tance; in three dimensions, the 2d surface states consist 
of single helical Dirac cones, and gapping them leads to 
an anomalous surface Hall response and the so-called ax- 
ion magneto-electric response 11 . The edge and surface 
states are at the heart of any possible technological ap- 
plication for topological insulators, including proposals 
to use them as a platform for Majorana states and quan- 
tum computation 14-16 . 

Alongside ideal TIs, recently, the study of the effects 
of imperfections - which plague essentially all realiza- 
tions so far - has also picked up. Experimentally, al- 
ready the 2d heterostructure realizations show the lack of 
quantized transport for samples of intermediate length 4 . 
Furthermore, 3d topological materials are often doped 
away from the bandgap, exhibiting a rather large bulk 
conductance 17 ' 18 . Theoretical studies of doping showed 
that the metal obtained when the chemical potential of 



the TI lies outside the bandgap is characterized by a fi- 
nite, but non-quantized topological index 19-23 . Studies 
of disorder led to additional surprises: it was shown that 
disorder can induce topological behavior in some trivial 
insulators 24-26 , and (more importantly for the current 
work) that disorder may close the gap in topological in- 
sulators, leading to a metallic phase 27-30 . 

In our work, we focus on the disorder-induced metallic 
phase in a simple test case of the Kane-Mele-Haldane 
honeycomb model 1,31 . Naively, once the disorder de- 
stroys the gap of a topological insulator, a simple 
metal emerges. By first constructing the self-consistent 
disorder-averaged Green function, we investigate the 
disorder-averaged topological index as well as the fate 
of the edge states in the metallic phase. Our results in- 
dicate that, on both counts, the disorder-induced metal 
retains some signatures of its topological origin. Just 
like a doped TI, it has a non-quantized, but finite topo- 
logical index. Correspondingly, we show that the edge 
states also survive as resonances of the disorder-averaged 
Green function. In particular, the edges states main- 
tain their helical nature over a finite lifetime/mean free 
path, and then get absorbed into the bulk. Surprisingly, 
their velocity is strongly renormalized by the disorder, 
and vanishes at the transition between the insulating and 
metallic phase. 

As pointed out by several groups, next-nearest neigh- 
bor interactions in various lattices may induce a topolog- 
ical phase 32-36 , namely a so-called Mott topological insu- 
lator. Using our formalism, we sclf-consistcntly analyze 
the formation of such a phase in the honeycomb lattice. 
We show that, in close analogy to superconductivity in 
the presence of magnetic impurities 37-39 , disorder sepa- 
rates the appearance of topological order into two tran- 
sitions: one associated with the appearance of a broken 
symmetry and an order parameter, and one (at larger 
interactions or weaker disorder) associated with the ap- 
pearance of a spectral gap. The phase with a finite or- 
der parameter but no gap coincides with the topological 
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metal wc discussed before. 

Below we start by developing the self-consistent 
disorder-averaged formalism (Sec. II), and use it to find 
the critical disorder strength for closing the TPs band gap 
(Sec. III). We then proceed to analyze the topological in- 
dex, following Volovik's prescription (Sec. IV). Perhaps 
the most physical consequences of our discussion regard 
the edge states. In Sec. V wc derive the properties of the 
edge states, assuming they correspond to a pole or reso- 
nance of the bulk disorder-averaged Green function when 
an edge is introduced. Before concluding, we discuss the 
possibility of an interaction-induced gapless topological 
phase in the disordered honeycomb lattice, and identify 
the parameter regime where this occurs (Sec. VI). 



where D is a high-energy cut-off of order of the band- 
width. 

One may distinguish two types of disorder, 
(o-icr z (7i)(sjS z Sj)(TkT z Tk) = ±a z s z r z . The strongest 
effect on the gap is due to impurities that obey 
((Ticr z (Ti)(sjS z Sj)(TkT z T k ) = (t z s z t z which is the case, e.g., 
for relatively smooth potential disorder, i = j = k = 0. 
Therefore, in the following, wc will restrict our attention 
to this case. 

To obtain the effective disorder-averaged Green func- 
tion for this problem, we note that the self-energy renor- 
malizes uj — > uj and A — > A. The self-consistent Born 
approximation is obtained by substituting the renormal- 
ized Cj and A in Eq. (7). The full Green function, 



G~ 



E, can then be written in the form 



II. DISORDER-AVERAGED GREEN 
FUNCTION 

We consider two copies of the Haldane Hamiltonian 
for a single-spin anomalous Hall state 31 or, equivalently, 
the Kane-Mele model for spin-1/2 fermions 1 . As we are 
interested in the low-energy physics, we restrict ourselves 
to the Hamiltonian near its Dirac nodes, 



% = vpk x a x T z — vpkyCTy — Aa z s z T z . 



(1) 



Here <j a are Pauli matrices acting in pseudo-spin (sublat- 
tice) space, while r 2 and s z are Pauli matrices operating 
in valley and spin space, respectively. 

Our goal is to understand the properties of this model 
in the presence of disorder which, if sufficiently strong, 
destroys its gap. Impurity scattering is described by 
adding 



(2) 



with (Uijk(v)} = and 

(Uijk(r)Ui mn (r')) = jijk S(r - Tc')5 a 5 jm 5 kn , (3) 

to the Hamiltonian (1). Here (...) denotes disorder av- 
eraging. 

Within the Born approximation, the self-energy due to 
impurity scattering reads 

£(w) = ^2 lijk OiSjTu I (dk) G (uj, k) (TtSjTk, (4) 



where 



G (uj, k) = (ioj - v F (k x a x r z - k y a y ) + Aa z s z r z ) 1 (5) 

is the unperturbed imaginary-time Green function of the 
system and (dk) = dk x dk y / (2ir) 2 . Upon momentum in- 
tegration, one obtains 
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G(uj, k) = [iu - v F (k x a x T z - k y a y ) + Aa z s z T z J , (7) 
where 

(8a) 



uj = uj(1 — C In —) : 

D 2 , 



A = A(l + Chi — 
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(8b) 



with C, = 7ooo/(47rwp'), using the definition of jijk in Eq. 
(2). 

In the absence of a topological insulator gap, A = 0, 
the above analysis readily yields the mean-free time. At 
A = A = 0, analytic continuation of Eq. (8a) to real 
time yields 



u> = uj I 1 — ( In 
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(iuj)< 



(9) 



For uj = 0, we then obtain the mean free time 40 at the 
Dirac point, r = i/uj(uj = 0). Namely, 



r" 1 = De~zc. 



(10) 



III. GAPLESS TOPOLOGICAL PHASE 



As noted by several authors 22 ' 23,27-30 , disorder pro- 
duces a phase transition between the topological insu- 
lator phase and a gapless metallic phase. This can be 
seen most directly using the self-consistent Green func- 
tion, given by Eqs. (7) and (8). In fact, the problem of 
finding the renormalized gap of the TI is equivalent to 
finding the renormalized gap of a superconductor with 
pair-breaking disorder. As explained in a review paper 
by Maki 39 , upon analytic continuation to real frequen- 
cies, the renormalized gap may be identified as the largest 
uj which permits a real solution for the renormalized uj 
and A. This follows from the expression for the density 
of states, 



= -T-Im 
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(dk) Tr[G(w,k)] 
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which upon using Eqs. (7) and (8a) yields 

1 



v(w) 



(12) 



Performing the analytic continuation and rcscaling all 
energies by A, Eqs. (8) simplify to 



UJ 

A= V 



Cm 



A 2 



5- 1 = 1 + Cln 



8 2 -v 2 
A 2 



8 2 -v 2 



(13a) 
(13b) 



where <5 = A/A, v = w/A, and A = D / A. Combining 
the two equations, we can eliminate the logarithm to find 



A 



(14) 



Solving Eq. (13b) for v and inserting the solution into 
Eq. (14), we then obtain 



'S 2 - A 2 exp 



-ifi-i 



which restricts the values of 5 yielding uj G 
one hand, 



S > 5 n 



(15) 



On the 



(16) 



and, on the other hand, in order for the square root to 
be real, 



S < S max = A exp 



1 

2C 



1 



1 



(17) 



The gap is then a function of £ and given by the maxi- 
mum of uj for the range <5 m i n < S < S max . 

At a critical value of the disorder parameter £, the gap 
vanishes, and an insulator-metal transition occurs. This 
happens when <5 max = S m - m = 1/2 and the two zeros of uj 
as a function of 5 merge. One obtains 



— = Ae 2t '§ a p 
2 

or, equivalently, 



(2-1) 
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with r given by Eq. (10), i.e., r ■ I = Dexp[-1/(2C, 



gap; 



IV. TOPOLOGICAL INDEX 

We are interested in determining whether the system, 
though gapless, retains some topological properties. Let 
us start by considering the disorder-averaged topological 
index of the system. There are several definitions for the 
topological index in 2d 6,10 ~ 12 ' 41 . For our purposes, the 
most useful approach is the one of Volovik 10,11,41 which 
calculates the Berry phase directly from the Green func- 
tion. The expression for the topological index, see e.g. 
Eq. (48) in Ref. 41, for one spin species is given by 



C=^JduJ (dk) e ij Tr [Gd u G~ 1 Gd^ G~~ 1 Gd^ G~ 1 ] CT T 



(20) 



In order to obtain the disorder-averaged topological in- 
dex, we simply substitute the disorder-averaged Green 
functions 42 . Using Eq. (7), one obtains 

G^G- 1 = G (id^ + d u Ka zSz T z ^ , (21) 
e ij Gd ki G- 1 Gd k .G- 1 = iv 2 F Go y {<j z t Zi G} o y . (22) 

Evaluating the fc-integrals and the trace, the expression 
can be reduced to 



where 



G = -s z [ duj (udujA - AduU)) 1 ~— . 
tt J V / w 2 + A 2 



Thus, one obtains 



2 1 
G = s z arctan , 

7f ^min 



(23) 



(24) 



Cj(uj = 0) 



Re 



A(w = 0) 
or, using the mean free time, 

1 - 



D 2 1 



(25) 




(26) 



In the gapped region, r _1 < T~Jp = A/2, one obtains 
^min = and therefore the topological constant is G = 
— s z - By contrast, if A = 0, one obtains x m - m = 00 
and therefore the topological constant vanishes, G = 0. 
The gapless region at r" 1 > t^L, however, emerges as 
the most interesting. In this region G changes smoothly 
from —s z to (see Fig. 1). Close to the closing of the 
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gap, < t 1 - r gap < r gap , wc find 

5C = C + s z oc ^rgap/r-1, (27) 

whereas deep in the gapless phase, r _1 3> r^p we find 
that the index falls off as 

Ccxr/r gap . (28) 

C(T) 




' 1 2 3 4 5 

T 

FIG. 1. The topological index for s z = — 1 as a function of 
the scattering rate 1/r. As soon as 1/r exceeds 1/T gap = 
A/2, the topological index plunges from a quantized value 
to a nonquantized and monotonically decreasing, disorder- 
dependent value. 

A non-quantized topological index requires a bit of con- 
sideration. For a given realization, by definition, Eq. (20) 
yields a quantized value. However, in the gapless phase, 
this value would fluctuate between and — s z , depend- 
ing on the specific disorder realization. Upon averag- 
ing, one may, thus, obtain a non-quantized value. It has 
already been pointed out by Loring and Hastings 22,23 , 
that the average topological index of doped 2d TIs in- 
terpolates continuously between the respective quantized 
values in the topologically non-trivial and trivial phases 
as the chemical potential is varied and the valence band 
becomes partially filled. Non-quantized values have also 
been obtained in clean 3d topological insulators, when 
the chemical potential is within the valence or conduc- 
tion band 19 ' 20 . 

The calculation we presented above fits well with the 
idea of a partial response of a metal. It shows that the 
disorder-induced gapless state should be quite similar in 
its behavior to a doped topological insulator, with the 
chemical potential lying within one of the bands. The 
index itself, however, is not a measurable quantity. Be- 
low, we try to understand the physical ramifications of 
this partial topological behavior in terms of the edge state 
physics. 

V. THE FATE OF THE EDGE STATE 

In the search for topological properties of the disor- 
dered metal, several works direct us to the edge states. 



Recently Ref. 13 confirmed rigorously that a spatial jump 
of the topological index (defined by Eq. (20)) results in 
an edge state, i.e., a pole 43 of the Green function at the 
location of the jump. It is natural to ask what a non- 
quantized jump implies for such edge states? 

Indeed, given the connections above, we speculate that 
the non-quantized topological index implies an edge state 
resonance of the Green function, but with a finite lifetime 
due to hybridization with the continuum. 

To explore the edge physics, we start with the inverse 
of the disorder-averaged Green function Eq. (7) which 
may be associated with an effective Schrodinger equation. 
For the four-component wave function ip Sz of a given spin 
species, the effective Schrodinger equation reads 

(v F k x a x T z ~ v F k y a y - Ka z s z T^j ip Sz = dnp Sz . (29) 

Consider now that the system occupies the half-plane 
x < 0. To satisfy the boundary conditions, we have to 
find solutions of (29) that vanish on the edge, x = 0. We 
recall that the degree of freedom t z = ±1 indicates the 
valley, i.e., the momentum around which k x , k y in the 
Hamiltonian (29) are measured. The true momenta are 
Px = k x + T z Ko and p y = k y , where Kq = 47r/(3-\/3a) and 
a is the lattice spacing. For the edge state to obey the 
boundary conditions, it is necessary to mix wave func- 
tions from the two valleys. Namely, the generic form of 
the edge-state wave function, assuming it is given by a 
simple superposition of momentum states (albeit with 
complex momenta), is ip — ip Tz (x, y)e' lTzKox . 

This implies that, for the wave function to vanish at 
the edge, both components t z = ±1 of the wave func- 
tion must be described by identical sublattice (cr-space) 
spinors. 

To find such solutions, we separate the Schrodinger 
equation into its parts dependent and independent of r 2 , 
and require that both vanish. Namely, 

(v F k x a x - Aa z s z ^\ 4> Sz = 0, (30a) 
(u) + vpkyGy) <p Sz — 0, (30b) 

where <p s . is a two-component wave function in sublat- 
tice space. In order for these equations to be solvable 
simultaneously, k x and k y must obey 

k x = ±is z — and k y = =)=——, (31) 

vf v f 

while the corresponding wave functions take the form 

4> Sz = (i,±t) T /V2. 

The complex valued momentum k x describes the decay 
of the edge state away from the edge. The decay length 
scale is, therefore, given by L± = up /Re [A]. The sign 
has to be chosen such that Im[fc x ] < 0. As can be seen 
in Fig. 2, the state remains well localized on the edge. 
Namely, the transverse localization length only increases 
from L± — fp/A at £ = to L± = 2up/A in the gapless 
phase. 
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FIG. 2. Penetration of the edge state into the bulk: The trans- 
verse localization length of the edge state, L± (r) = w_p/Re[A] , 
in units of the bare decay length, L±(0) = Vp/A, is plotted 
as a function of disorder. It varies by a factor of 2. 



Vf Re/Im^] 



(b) 



The momentum k y describes the propagation along 
the edge. The appropriate choice of sign for k x yields 
k y = —s z Q/vp. Thus, as expected, the sign of k y and 
therefore the propagation direction is tied to the spin. 
To obtain the dispersion of the edge states, as well as 
their lifetime, we need to find the relation between k y 
and the real energy u>. We use Eqs. (8) to write 



u = u) [ 1 — C In — 



D 2 (2Cj-uf 



w 2 (A 2 - {2Co-lo) 2 ) 



(32) 



Since the edge state is helical, an imaginary part of the 
wave number k v can be interpreted as a finite lifetime 
of the state. To illustrate this, consider the simplest 
Schrodinger equation for a decaying helical mode: ^ — 



z ox 



di 



This is cquivalcntly solved either by tp 



cxp[iuj(t — x/v) — t/r] or ip = cxp[iu(t — x/v) — xj (vr)]. 
Which of these two solutions should be used depends on 
the specific situation, and, for a Green-function descrip- 
tion, the two should be equivalent. 

When the disorder is sufficiently weak, ( < Cgap, the 
dispersion decribed by (32) can be separated into two 
regimes. At energies below the renormalized gap edge, 
u> < -Egap < A, we find edge state solutions that do not 
decay, i.e., Im[w] = 0. By contrast, at energies co > E g&p , 
we find edge states with a finite lifetime, Im[a)] =/= 0. In- 
deed, the finite lifetime of the edge state is a consequence 
of the hybridization with the continuum which permits 
for the state to decay to extended bulk states. The dis- 
persion and lifetime of the weak-disorder edge states is 
shown in Fig. 3a. 

As the disorder increases, C > Cgap, the bandgap van- 
ishes. Now the entire edge state branch is hampered by a 
finite lifetime and propagation length, see Fig. 3b. The 
decay length in the propagation direction L|| is nearly 
independent of energy, but depends strongly on the dis- 
order strength. By setting ui = Rc[w] = 0, one can obtain 
L|| = vp /\m[Cj\. We find, relying on the identification of 
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FIG. 3. Dispersion of the edge states: k\\ = Re[d)]/wir (full 
line) and inverse decay length L^ 1 = Im[o;]/n.F (dashed line) 
versus oj. (a) Gapped phase (C, = 0.8 Cgap)- (b) Gapless phase 

(C= 1.2 Cgap). 



k y in Eq. (31), 




(33) 



where r is the mean free time as defined in Eq. (10). The 
result is shown in Fig. 4. At ( = Cgap, the decay length 
diverges, L|| — > oo. 

Quite interestingly, the survival length of the edge 
state is closely related to the topological index calculated 
in Sec. IV. Namely, 



2 in A 

C = s z arctan — - — . 

7T 2vp 



(34) 



This gives more ground to the speculation that a non- 
quantized topological index will always be associated 
with a finite lifetime for the edge states. 

Another interesting property of the edge states is the 
propagation velocity, v — vp{doj / du). In the limit u — > 
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FIG. 4. Inverse decay length along the edge, L,, 1 , in units of 

the inverse bulk mean free path, i _1 = 1/(vft), as a function 
of disorder. L7 1 is zero in the gapped phase and approaches 



V(T) 



(a) 




r(r)r 



1.0 - 



0, straightforward manipulation of Eq. (32) leads to 



v 

V F 



1 - 2(ln 



D(2- 



2vi,T-'L: 



- 



C<Cc, 

OCc 



(35) 



The result is depicted in Fig. 5a. Approaching the tran- 
sition from both sides, we find that the velocity vanishes 
when C = Cgap- In the gapped regime, we identify the 
transcendental equation defining v for C < Cgap as reduc- 
ing to the gaplessness condition (18) upon setting v = 0. 
Indeed, the edge velocity is suppressed throughout the 
weak-disorder regime at energies tending to the gap from 



below, 



E gap . In the gapless regime, C > Cg 



the 



denominator 44 in Eq. (35) is regular at £ = Cgap- Thus, 
since L^ 1 vanishes at the transition, so does v. 

With the velocity and the decay length, we may now 
identify the decay rate of the edge states in the gapless 
regime as a function of disorder, 



T = — = 



2v 3 f t 2 



1 



L 



c- 



(!*■) 



(36) 



Thus the lifetime ~ r -1 diverges at the point where 
the gap closes, and follows the third power of the decay 
length as disorder further grows (the result is depicted in 
Fig. 5b). 



VI. TOPOLOGICAL MOTT METAL 

Topological behavior may also emerge due to electron- 
electron interactions rather than a native spin-orbit cou- 
pling and band inversion 32-36 . Under these circum- 
stances, it is particularly interesting to ask whether a 
gapless topological phase can exist when the disorder is 




FIG. 5. Properties of the edge states as a function of the 
scattering rate r _1 . (a) Effective velocity v of the edge states 
in units of the Fermi velocity, (b) Decay rate F of the edge 
states in units of the scattering rate r -1 . As the velocity 
vanishes at Cgap, F increases at a much slower rate than L^ 1 
in the gapless phase, namely F oc 3 . 



strong. The situation is rather reminiscent of the compe- 
tition between magnetic impurities and s-wave supercon- 
ductivity. It is well known that magnetic impurities may 
induce a phase which is s-wave superconducting, albeit 
gapless 37-39 . Namely, the disorder suppresses both the 
order parameter and the spectral gap, but the two do 
not vanish simultaneously. In the gapless phase, while 
the material has a superconducting stiffness, its critical 
current is suppressed. In this section, we will establish 
that interacting systems can host a gapless topological 
phase as well, albeit in a rather narrow range of interac- 
tions and disorder. 

The idea of an interaction-induced topological insula- 
tor, i.e., a topological Mott insulator, was first explored 
in Rcf. 32 for the case of a honeycomb lattice. It was 
shown that a next-nearest neighbor interaction may open 
a gap which has opposite sign at the two Dirac nodes of 
the band structure. Using the Hartrce-Fock approach 
this can be analyzed self-consistently, which is the ap- 
proach wc will pursue below. Later work 33 demonstrated 
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that various lattice structures may also be susceptible to 
competing symmetry-broken phases, most notably the 
Kekule phase in the honeycomb lattice. We will stir clear 
of the possibility of competing phases, and concentrate 
on the range of parameters most conducive to the forma- 
tion of topological order. 

In particular, we will consider a honeycomb lattice in 
2d with next-nearest neighbor repulsion of strength V. 
Our starting point is the tight-binding Hamiltonian, 



H = -t 



(37) 



+^E fecU-l) fee] 



^3° - 1 



To this Hamiltonian, we will add disorder as explained 
in Eq. (2). 



A. Self-consistent Hartree-Fock analysis 



nearest neighbor bonds. The next-nearest neighbor hop- 
ping may acquire an imaginary expectation value which 
will then result in the Haldanc model for each spin fla- 
vor. A next-nearest neighbor interaction naturally leads 
to such a broken symmetry 32 ' 33 . The interaction term 
may be decoupled using a Hubbard-Stratonovich trans- 
formation, yielding the mean-field Hamiltonian, 



H = - t E ( c Uv+n.c.)-2y Xtfclc ja >, (38) 



together with the self-consistency equation, 



Xij ( c jV Ci < 



(39) 



The symmetry breaking that leads to topological be- 
havior in the honeycomb lattice occurs on the next- 



J 



After Fourier transformation, the Hamiltonian in the 
sublattice (pseudo-spin) basis takes the form 



H 



E 

k.k'icr,^' 



a 



krr 



A CTCr 



(k,k') 

^k.k^, 



k,k'<W 



-ika. 



! <5k.k'<5, 



k,k'CW' 



AJT (k,k') 



(40) 



where A^ = IVXx^j and X = A, B is the sublattice 
index. Furthermore, a; (with I = 1,2,3) are the lattice 
vectors of the honeycomb lattice, connecting the A and 
B sublattices. 

At low energies, we assume that the interaction opens a 
topological insulator gap without breaking time-reversal 
symmetry, and therefore A has the form 



A° x ° (k,k') = Asign^X-l)^, 



k.k' 



(41) 



with {-l) A / B = ±1. Note that s\gn(k x ) inherits the 
valley index, and therefore in a continuum picture will 
be replaced by r 2 . Indeed, the low-energy Hamiltonian 
of the system, upon linearization of the dispersion near 
the nodes, reads 



H = v F k x a x T z - v F k v <7 v - A<t z s z t z 



(42) 



which is the same as Eq. (1). However, this Hamiltonian 
is now supplemented by the self-consistency condition for 
A. Assuming that the low-energy part of the dispersion 
is dominant in determining the self-consistent bond order 
parameter, the self-consistency condition takes the form 



A = V J ^ J (dk) Tr [G(u, k)a zSzTz ] (43) 
with G(w,k) given by Eq. (7). 



B. Finding the critical disorder and interaction 

Disorder now has two effects. As before, it suppresses 
the spectral gap, but, at the same time, it also affects the 
order parameter: for any interaction strength which leads 
to the topological phase, adding sufficiently strong disor- 
der will destroy that phase. The condition for the spec- 
tral gap to vanish remains the same as before. Namely, 



Cgap — 



21n(2f) 



(44) 



see Eq. (18). However, A = A(C gap ) has now to be deter- 
mined self-consistently. Substituting the Green function, 
Eq. (7), into the self-consistency condition (43), we find 



A = 8V 



2^ 



{dk) 



A 



uj 2 + v 2 F k 2 + A 2 



(45) 



which, upon performing the momentum integration and 
using Eqs. (8), reduces to 



c 

CA 



du (1 -6), 



(46) 



with c = 2VD/{Tr 2 Vp), and A = D/A as defined earlier. 
Note that 5 and uj are related through Eq. (15). 
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Fig. 6 shows the order parameter Ao in the absence of 
disorder. To obtain the critical disorder strength, £ c , for 
which the topological phase is destroyed, we replace the 
w-integral with an integral over S and let 6/ A — > 0. The 
result is a transcendental equation that connects c and 

Co, 

1/2 

The critical scattering rate 1 /r c as a function of c is plot- 
ted in Fig. 6. 




FIG. 6. Order parameter in the absence of disorder, Aq/D, 
and critical scattering rates, 1/(Dt c ) (vanishing of the order 
parameter) and l/(Dr ga p) (vanishing of the spectral gap), as 
a function of interaction strength, c. 



C. Disorder threshold for gaplessness 

To obtain the disorder strength, Cgap, where the spec- 
tral gap vanishes, we need to combine Eqs. (43) and (46). 
We find another transcendental equation for Cgap; 

!_ 1 

1 = J dS V / e -C^ 1 p(5- 1 -2)_ 4(5 2. (48) 

Sap 1/2 

The result is plotted in Fig. 6. Furthermore, the width 
of the gapless regime is shown in Fig. 7. The larger 
c and thus Ao, the larger the regime where the order 
parameter survives while the gap is closed. In the gapless 
regime £ gap < £ < £ c , the above considerations about the 
topological index and edge states remain valid, i.e., the 
system is a topological Mott metal. 

VII. CONCLUSIONS 

We have studied the disorder-induced gapless phase 
of a 2d topological insulator system. In particular, we 




FIG. 7. Region of gaplessness at T — 0: (t c _1 — r^^/r" 1 as 
a function of 1/(Dt c ). 



have concentrated on (i) the Kanc-Mcle model, which is 
intrinsically topological, as well as (ii) an example of a 
Mott topological insulator, namely graphene with next- 
nearest neighbor interactions that induce a topological 
phase. 

Our main conclusion is that the disorder-induced gap- 
less system remains topologically non-trivial. Naively, 
one would think that once the gap of a topological in- 
sulator is closed, its topological properties also vanish. 
Nevertheless, it appears that the gapless phase retains ev- 
idence of its topological origin. Namely, the gapless phase 
is characterized by a finite, yet non-quantized disorder- 
averaged topological index. This is reminiscent of doped 
topological insulators for which Refs. 22 and 23 showed 
that the topological index changes continuously between 
1 and when the chemical potential is sweeped between 
the top and bottom of the valence band. 

In addition, we find that the gapless phase still sup- 
ports helical edge states, though with a finite lifetime and 
a renormalized velocity. This is consistent with the find- 
ings of Ref. 21, where a bulk parasitic metal overlapping 
with edge states was considered in the clean case. There 
the individual orbitals of the edges states which overlap 
in both momentum and energy with bulk states arc split 
into two states, with energies above and below the para- 
sitic metallic band. However, these states leave a 'ghost' 
behind: an edge state with a finite lifetime, which ap- 
pears as a resonance in the Green function of the system, 
although not being an exact eigenstate of the Hamilto- 
nian. Furthermore, the finding of a finite lifetime limiting 
the manifestations of topological effects may be related 
to the 'parity diffusion mode' found in Ref. 27 (see, e.g., 
Eqs. (4) and (66) ibid., where the IR cutoff indicates a 
finite lifetime) . Since that work considered disordered 3d 
TIs, a direct comparison, however, is not possible. 

The properties of the edge states/resonances have ad- 
ditional features worth mentioning. First, it appears that 
a close link exists between the decay length of the edge 
resonances and the non-quantized topological index, as 
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shown in Eq. (34). Next, disorder strongly affects the 
edge state velocity. Weak disorder suppresses the edge 
velocity, until it vanishes at the disorder strength where 
the gap closes. This is essentially dictated by the fact 
that the helical edge branch is confined between the top 
of the valence band and the bottom of the conduction 
band. Interestingly, the velocity picks up as a function 
of disorder in the gapless regime, as shown in Fig. 5b. 
The suppression of the velocity when £ ~ £ gap implies 
that an electron tunneling into the surface will remain 
near the point of tunneling for an extended time. This 
effect could be observed through an enhanced zero-bias 
anomaly due to a temporary Coulomb blockade 45-47 . 

If the topological phase is induced by interactions, as 
considered in the last part of the paper, the disorder af- 
fects both the excitation gap and the order parameter. In 



analogy with gapless superconductivity, we have shown 
that there is a range of disorder strengths where a gapless 
topological phase is realized. In this regime, the disorder 
fully suppresses the gap, whereas a finite order param- 
eter prevails up to a critical disorder strength. Within 
this range, the material will have the same properties as 
the non-interacting disorder-induced gapless topological 
phase. 
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